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Abstract: This paper presents the application of a frequency-domain reflectometry (FDR) 
sensor designed for soil salinity assessment of sandy mineral soils in a wide range of soil 
moisture and bulk electrical conductivity, through the determination of soil complex 
dielectric permittivity spectra in the frequency range 10-500 MHz. The real part of 
dielectric permittivity was assessed from the 380^40 MHz, while the bulk electrical 
conductivity was calculated from the 165-325 MHz range. The FDR technique allows 
determination of bulk electrical conductivity from the imaginary part of the complex 
dielectric permittivity, without disregarding the dielectric losses. The soil salinity status 
was determined using the salinity index, defined as a partial derivative of the soil bulk 
electrical conductivity with respect to the real part of the soil complex dielectric 
permittivity. The salinity index method enables determining the soil water electrical 
conductivity value. For the five sandy mineral soils that have been tested, the relationship 
between bulk electrical conductivity and the real part of dielectric permittivity is essentially 
linear. As a result, the salinity index method applied for FDR measurements may be 
adapted to field use after examination of loam and clayey soils. 
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1. Introduction 

Accurate and reliable estimation of soil salinity, defined as the electrical conductivity of soil water 
extract or saturated soil water extract, is a very important issue, especially in arid regions, where 
salinity of the soil may increase significantly and pose a danger to plants [1,2]. Soil salinity 
measurement is rather complicated as it is often affected by factors producing measurement errors, 
e.g., temperature, moisture and texture of the soil. Therefore, the development of new measurement 
methods and equipment for evaluating soil salinity status is still an active field of research. The ideal 
measurement tool for soil salinity status should read the majority of influencing components at the 
same time and in the same location as quickly as possible to register the momentary values of interest. 
This is done by incorporating various sensors in a single unit for insertion into the measured 
material [3,4] or by selective analysis of sensors' outputs for discrimination of various quantities like 
soil water content and electrical conductivity by time domain reflectometry [5]. The discrimination 
techniques sometimes require increased accuracy of measurement and additional calibrations, like 
calibration of TDR probes for low soil water content measurement [6]. Application of indirect 
measurement sensors based on soil dielectric properties, i.e., TDR and FDR sensors, has opened new 
possibilities through the ease of integrating the sensing elements of soil water content, electrical 
conductivity and temperature in one small measurement probe. The concept of the soil salinity 
index X s , describing soil salinity status, was introduced by Malicki and Walczak [7] as one of many 
applications of the TDR technique for nondestructive and simultaneous measurement of water state 
and its transport in soil [8]. The salinity index was defined as a partial derivative of the soil 
bulk electrical conductivity C b with respect to the soil bulk or apparent dielectric permittivity s b 
(Equation (1)), where both variables were determined from the TDR waveform at the same time and 
on the same soil volume: 

dC b 

Soil bulk dielectric permittivity and electrical conductivity were determined from the velocity of the 
TDR pulse and its attenuation, respectively, when traveling along the TDR parallel waveguide inserted 
in the soil [5,9]. Malicki and Walczak [7] showed that the X s value depends primarily on soil salinity 
and soil texture. Soil salinity can thus be uniquely determined from the salinity index X s and the sand 
content of the soil. Their experimental evidence showed that the soil salinity index as a function of soil 
electrical conductivity of the electrolyte is independent of the soil volumetric water content for values 
of water content above 0.2. 

The application of the TDR technique in calculating the salinity index includes simplifications that 
can affect the final result. The value of s b approximates the real part e' of the complex dielectric 
permittivity and is dependent on the frequency / of the applied external electric field. Also, the 
electrical conductivity, C b , calculated from the attenuation of the TDR pulse, does not take into 
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account the dielectric loss e d . The complex dielectric permittivity e*, presented by Equation (2), 
comprises the real part that describes the energy storage in dielectric medium, and the imaginary part 
that measures the energy losses. These losses may be caused by dissipation of energy during 
polarization processes (dielectric loss) and by electrical conductivity of the material: 



where e 0 = 8.85 x 10~ 12 F m _1 stands for dielectric permittivity of free space. 

The dependence of the complex dielectric permittivity of soils on the frequency of the applied 
electrical field has been discussed in literature; however, the applied frequency was not defined [10], 
the interpretation was limited to 50 or 150 MHz [11,12] or the frequency spectrum was obtained from 
the TDR waveform [12]. 

It was shown [10] that the separation of real and imaginary parts of the complex dielectric 
permittivity would remove the influence of the conductivity effect on the water content calibration of 
the FDR technique compared with TDR. The application of the broadband FDR technique allows 
measuring the real and imaginary parts of the complex dielectric permittivity of soil separately and 
independently. Consequently, for water content measurements it is possible replace the bulk dielectric 
permittivity (the product of TDR) with the real part of the dielectric permittivity determined by FDR in 
the frequency range near 400-500 MHz, where the influence of electrical conductivity on the real part 
of dielectric permittivity is negligible. Another reason for developing FDR devices for measuring soil 
moisture and salinity is that the TDR meters require the application of steep pulses with very short 
rising times (at least several times less than a nanosecond), which significantly contributes to the cost 
of these devices. 

The aim of the paper is to determine the bulk electrical conductivity and soil pore water electrical 
conductivity of selected soil samples based on independent measurements of the real and imaginary 
parts of the complex dielectric permittivity using the FDR technique in the frequency range of 10 MHz 
to 500 MHz. The received FDR-based linear salinity index is compared with the Malicki and 
Walczak [7] TDR-based salinity index presented in Equation (1). The quadratic FDR-based salinity 
index model, which takes into account possible quadratic contributions to the C b vs. s' relation, is 
presented and its performance is compared to the linear model. The possibility of obtaining the soil 
pore water electrical conductivity value from a single FDR measurement of C b and s', necessary for 
adapting the FDR-based salinity index approach to field use, is also evaluated. 

2. Materials and Methods 

The tested material included five mineral soils that were air dried and put through a 2 mm sieve. 
The basic physical characteristics of these materials are given in Table 1 . The soils were collected from 
the Ap layer, 20-30 cm below the surface. The material dry bulk density was determined in the 
laboratory in 10-fold repetitions by weighing soil samples of known volumes, drying them at 105 °C 
for 24 hours and again weighing [13]. The respective values from Table 1 are calculated for a 
confidence level of 95%. Having determined the material bulk density, it was possible to calculate its 
mass in a measurement container of 120 cm 3 volume (cylinder; diameter: 4.6 cm, height: 7 cm), which 




(2) 
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helped to prepare appropriate soil samples with variable water content ranging from approximately 
10% to near saturation. The soil material of various moisture levels was packed into the containers in 
small portions and pressed with a 0.2 kg rubber hammer to achieve homogenous density distribution in 
the soil sample [7]. There was no intention of achieving predefined values of soil sample density, as 
the experiment tested only the idea of the FDR-determined salinity index. Five soil materials were 
tested. Each soil was wetted with four solutions (distilled water and three KC1 solutions, presented in 
Table 2), to seven soil moisture content levels. There were 35 pots with soil samples for each wetting 
solution, giving a total number of 140 pots. FDR measurements were performed three times on each 
sample. 



Table 1. Selected physical characteristics of tested material; RSD stands for relative 
standard deviation (standard deviation/mean value). 



Name 


Bulk density 


Specific surface a 


Saturation water content by mass (%) 


Texture (FAO) 


" (%) 


(kg-m 3 ) 


(mV) 


Average 


RSD 


sand 


silt 


clay 


568 


1,415 ±3 


34 


37.2 


0.03 


58 


31 


11 


589 


1,697 ±4 


10 


22.4 


0.04 


88 


11 


1 


601 


1,380 ±4 


31 


38.7 


0.02 


60 


34 


6 


605 


1,635 ±7 


10 


23.8 


0.07 


95 


4 


1 


loess 


1,382 ±4 


30 


37.6 


0.01 


55 


29 


16 



Water vapor adsorption method [14]. 
b Data taken from Glinski etal. [15]. 



The values of saturation water content by mass for the selected soils were determined using 
capillary rise. Three samples of each soil of a known dry weight placed in measurement containers 
with holes at the bottom, were gathered in separate plastic bowls with porous bottom sides. The 
containers with soil, covered to minimize excessive evaporation, were immersed for 72 hours in 
distilled water allowing capillary rise to achieve saturation. Then they were weighed to determine the 
mass of water saturating each soil sample. Finally, after simple calculations it was possible to 
determine the mass of water for adding to the 120 cm volume of air-dry soil to achieve soil samples of 
a desired water content, from 10% to near saturation. 

On the basis of generally available conversion tables and pilot measurements of soil electrical 
conductivity performed by a TDR meter, three KC1 solutions and distilled water were prepared for 
wetting soil samples (Table 2). Application of KC1 solutions of various electrical conductivity C s 
enabled changing the real and imaginary parts of dielectric permittivity of the tested material to 
achieve the assumed range of variability of the soil electrical conductivity of 100-400 mS m~ , and of 
soil water content — from air dryness to near saturation. The samples of chosen soils fully filled 
120 cm 3 containers (plastic cylinder 4.6 cm in diameter and 7 cm in height) equipped with a sealed 
cover, and wetted by distilled water and the KC1 solutions. 

The material samples in containers were mixed with distilled water and KC1 solutions, covered with 
a sealing cap, weighed and conditioned at 40 °C for 72 hours to ensure uniform water content in the 
sample volume. Then, after leaving the soil containers for several hours at room temperature, the filled 
containers were again weighed to make weight corrections caused by possible water evaporation loss. 
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Next, the containers were opened to perform FDR measurements of the complex dielectric permittivity 
of the material. Each FDR measurement was made three times by inserting the probe rods in various 
locations of the material in the container. All measurements were made in a laboratory with a 
controlled temperature 21 ± 1°C. 

Table 2. Parameters of KC1 solutions applied for wetting the measured material (values at 
temperature 21 °C). Assumed C s values were calculated from molarity [16] and confirmed 
by measurements. 



Label of solution 


KC1 molarity (mol dm 3 ) 


Assumed C s (mS-m ') 


1 


0.000 




2 


0.042 


490 


3 


0.093 


1,080 


4 


0.150 


1,680 



2.1. FDR Probe 

The applied FDR sensor, measurement details and calibration techniques were described earlier by 
Skierucha and Wilczek [17]. The probe consisted of a parallel waveguide with two steel rods, similar 
to probes used with the TDR instruments developed at the Institute of Agrophysics PAS in Lublin, 
Poland (easytest.ipan.lublin.pl), but shorter, i.e., 3 cm against 10 cm. The probe with such short rods 
applied with a TDR meter would require even shorter rising time of the pulse, which is possible in 
extremely sophisticated and expensive TDR devices. With a probe of this type, the measured quantities 
pertained to a greater volume of the sample under test than in the case of an open-ended coax probe [18]. 
The soil is a highly inhomogeneous, multiphase material and testing of it provides many challenges, as 
the measurement result taken at a given point in the sample could differ from the result obtained at 
another point. Therefore, the open-ended coax probe is impractical in this case, as it would introduce 
significant errors due to the inhomogeneity of the material. A vector network analyzer (VNA), type 
ZVCE from Rohde and Schwarz, was used for measuring the complex reflection coefficient 5 1X of the 
signal reflected from the probe inserted into the sample. This reflection coefficient is defined as: 

_ Z p (e*)-Z c 

A " Z p {e*)+Z c V* 

where Z p is the impedance of the probe and Z c = 50 H. The impedance of the probe depends on the 
complex dielectric permittivity of the material in which the sensor was inserted. Therefore, the 
measured reflection coefficient is also a function of the dielectric permittivity of the sample. It is 
possible to obtain the inverse relation: s* = /(Sn), which allows determining the complex dielectric 
permittivity from the measured reflection coefficient for each applied frequency. For calibration 
purposes, measurements of short circuit, open (air) and acetone were used, which reduced the 
measurement errors [19]. 
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3. Results and Discussion 

3.1. Complex Dielectric Permittivity 

The FDR measurement technique allowed determining the real and imaginary parts of the dielectric 
permittivity directly and independently. Figure 1 presents the real and imaginary parts of the dielectric 
permittivity measured for a sample of soil no. 601 wetted with distilled water to approximately 50% of 
the saturation water content. The other soil samples were tested in a similar way. On the graph one can 
notice artifacts located near frequencies of 150, 205, 360 and 500 MHz for both parts of s*. They are 
related to resonances in the experimental set-up, due mostly to the length of the coax cable connecting 
the sensor with the VNA; the details and discussion are presented in [17]. This was later confirmed by 
introducing a magnetic shield on the coax cable, which substantially decreased the effect of the 
artifacts. The frequencies for which the artifacts occur will be excluded from further analysis, without 
negatively influencing the obtained results. 

Figure 1. Real and imaginary parts of the complex dielectric permittivity of a sample of 
soil no. 601 wetted with distilled water to approximately 50% of saturation water content. 

25 t 




o -I 1 1 1 1 1 1 

0 100 200 300 400 500 600 

f /MHz 

The real part of the dielectric permittivity is strongly related to the soil water content, namely the 
square root of the real part of the dielectric permittivity (the refractive index) depends linearly [20] on 
the volumetric water content, measured by the standard bulk density plus thermogravimetric soil water 
content method. This dependence for results obtained by the FDR probe used in this experiment is 
confirmed in [18]. To minimize the effect of electrical conductivity of the sample on the real part of 
dielectric permittivity, the values of the real part of the dielectric permittivity were taken for the 
frequency range 380^140 MHz [17]. 
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3.2. Electrical Conductivity of Soil Samples 

The imaginary part of the complex dielectric permittivity may be used to infer the bulk electrical 
conductivity of the sample. According to Equation (2), the imaginary part of the complex dielectric 
permittivity e* is related to bulk electrical conductivity of the sample C b 

[m(£ " ) = £ " + 2^7 (4> 

One can multiply both sides of the above equation by the frequency, so that the whole relation is a 
linear function of /: 

Im(£*) -f = e d f + 



2ns 0 



(5) 



Assuming that e d is not frequency dependent in the analyzed frequency range up to 500 MHz, 
one can find the values of C b by fitting a straight line into a plot of the above function, as shown in 
Figure 2. The regression equation, coefficient of determination R 2 and standard error of regression a, 
defined as the square root of the sum of the squared residuals divided by the number of degrees of 
freedom, were presented on the graph. To calculate bulk electrical conductivity with the artifacts 
removed, frequencies from 165 to 180 and from 245 to 325 MHz were selected. In this frequency 
range, which is narrow enough and distant from relaxation frequencies of various polarization 
mechanisms which may occur in the tested soils, the assumption that the dielectric loss e d does not 
depend on frequency is reasonable and gives good fits, as presented by the example in Figure 2. 

Figure 2. Determination of bulk electrical conductivity of a sample of soil no. 601 wetted 
with distilled water to 50% of saturation water content; C b is directly proportional to the 
intercept, according to Equation (5). The regression equation, R 2 and standard error of 
regression are given on the plot. The standard error of determination of C b is about 2%. 



950 -) 

900 - 

850 - 
800 

^J, 750 - 

E 

~ 700 H 

650 
600 - 
550 - 
500 




1.27x + 508.45 
= 0.99, 0=7.69 



100 



150 



200 250 
f /MHz 



300 



350 



This procedure was applied to determine the bulk electrical conductivity C b of all tested 
soil samples. 
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3.3. Salinity Index and Electrical Conductivity of Soil Water 

In order to test the concept of salinity index introduced by Malicki and Walczak [7] with respect to 
the FDR measurements of the real and imaginary parts of the dielectric permittivity of soils, the 
electrical conductivity C b of each sample versus the real part of dielectric permittivity s' = Re(£*) was 
plotted for all moistening solutions and a straight line was fitted through each set of data points. The 
results are presented on the left panel of Figure 3. The regression equation, R 2 and standard error of 
regression a, defined as in the previous section, were presented on the graphs. For all points, error bars 
are present. The vertical error bars represent the standard error of determination of C b . Taking into 
account these errors and high values of R 2 (as presented on the graph), it transpires that for all tested 
samples the relations between C b and e' are essentially linear. Thus, the salinity index defined as: 



is actually the slope of the fitted line in the C b vs. e' graph and, for linear C b vs. e' dependence, is 
independent of the water content. The difference between Equation (6) and Equation (1), which 
presents the salinity index defined for TDR measurements, is the use of the real part of the dielectric 
permittivity instead of the apparent dielectric permittivity. In addition, the electrical conductivity in 
Equation (6) was determined from the frequency spectrum of the imaginary part of the dielectric 
permittivity as described in Section 3.2, while in Equation (1) the quantity C b actually denoted bulk (or 
apparent) electric conductivity determined from attenuation of the TDR pulse with the assumption that 
the dielectric loss s d = 0. 

The relations between the salinity index and moistening solution conductivity C s for the examined 
soils are presented on the right panel of Figure 3. Each value of X s was determined from the slope of 
the regression equations of C b against s' , which were shown on the left panel of Figure 3. It was found 
that X s depends linearly on the electrical conductivity of the moistening solution C s . When the salinity 
index X s and the slope I of the X s vs. C s relation are known, one can calculate the conductivity of soil 
water from the following formula: 



One may notice that for the samples wetted with distilled water the salinity index is equal to some 
initial value X SI 0, due to some residue conductivity of ions dissolved from dry soil. As was shown 
in [7], this residue conductivity C r may be found by extrapolating the X s vs. C s relation to the 
horizontal axis. Then: 



Therefore, one may expect that the electrical conductivity of soil water C w is a sum of the residue 
electrical conductivity and the conductivity of the moistening solution: 




(7) 




(8) 



C r + C s 



(9) 



This assumption will be tested in a subsequent part of this paper. 
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Figure 3. On the left: bulk electrical conductivities of all measured soil samples vs. the real 
part of the dielectric permittivity. On the right: salinity index X s calculated as a slope of a 
linear C b vs. e' relation, with respect to the conductivity of the moistening solutions C s . 
Regression equations, R 2 and standard errors of regressions are given on the plots. 
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3.4. Salinity Index Method for Field Use 

To calculate the salinity index using the method presented above, it is necessary to take a series of 
measurements of the same soil moistened with the same solution to various water contents. Obviously, 
this procedure has little practical use, since a measure of soil salinity applicable for field conditions 
should provide an accurate estimate based on a single measurement of a single soil sample. However, 
similarly to the method shown in [7], if the relation C b vs. e' is linear, one need not take a series of 
measurements of samples of different water contents of the soil under question to determine the value 
of the salinity index — the measurement of complex dielectric permittivity of a single sample will 
suffice. It transpires that the lines from the left panel of Figure 3, fitted into the (V, C b ) data points, for 
a given soil and all applied solutions, cross at certain limiting values of the real part of the dielectric 
permittivity and bulk electrical conductivity. These limiting values are denoted s t and C h respectively. 
Once these values are known for a given soil, it is possible to calculate the salinity index by 
determining the value of the partial derivative from an appropriate difference quotient: 

_ dC b AC b _ C b - Cj 

Xs ~j^ = jf~v^t; (10) 

with the use of the values of e' and C b measured from a single sample of previously unknown salinity 
and water content. Therefore, combining Equations (7) and (10), the electrical conductivity of soil pore 
water may be calculated from the following relation: 

C b Ci 

To apply the formula above in the field, it is necessary to know the values of £/ and C t for a given 
soil. The values, obtained in this experiment for the five tested soils, depend on the soil properties. It 
was found that s { depends on the specific surface s of the soil under question, with regression equation 
given by £/ = 0.02s + 3.44 with R 2 = 0.99, where s is given in units from Table 1. On the other 
hand, for the tested soils C t is affected by the clay content c. The appropriate regression equation is 
C; = 0.65c + 3.23, with R 2 = 0.89 and clay content given in units from Table 1. However, before 
these equations can be used in the field for an unknown soil, it is necessary to perform laboratory 
measurements of a greater number of soils with different texture and other properties to obtain more 
reliable regression equations. Because of the limited number of the tested soils, the regression 
equations on s { and C/ are included only to present the variability of these parameters with the soil 
texture and are not used in further calculations. The accuracy of the values of £ { and C/ is especially 
important for soils which are either very dry, or have low bulk electrical conductivity. This is because 
from the form of the last term of Equation (10) it transpires that for samples with values of C b and e' 
close to C 1 and s h respectively, the error of the final value is the most affected by the errors of all these 
four quantities. 

To calculate the electrical conductivity of soil water C w from the salinity index X s in the field 
conditions from Equation (1 1), the value of the parameter I needs to be known beforehand too. As was 
the case with e { and C /5 the parameter I also depends on the soil properties. For the five tested soils, the 
best fit was obtained for a straight line of equation I = —0.00008c + 0.009, where c is the clay 
content, as for C h However, R 2 equals only 0.56. It is obvious that in order to minimize the error of 



Sensors 2012, 12 



10900 



C w determined in field conditions, additional soils of different properties (especially loam and clayey 
soils) should be tested in order to achieve better regression equations. 



Figure 4. Soil water conductivities C w versus e' = Re(£*) for linear and quadratic salinity 
index models for all moistening solutions. Solid lines represent C w values calculated with 
the use of X s obtained from the slopes of C b vs. e' linear relations. 
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3.5. Comparison of Linear Models of Salinity Index 

On the left panel of Figure 4, the values of soil water electrical conductivity C w for all tested soils 
are shown. The "linear model" phrase at the top of the figure means that the values of C w were 
calculated with the assumption that the relation C b vs. s' is linear as in the left panel of Figure 3. 
Possible quadratic effects are taken into account in the "quadratic model" elaborated in the next 
section. The straight lines on all graphs represent the values of C w (these values are also written on the 
right side of the lines) calculated from Equation (7) with the use of the salinity index X s obtained from 
the regression equations from the left panel of Figure 3 and I values determined from the slopes of the 
regression lines from the right panel of Figure 3. Since the salinity index calculated from the slope of a 
straight line fitted through the data points obtained from the whole series of samples moistened to 
various water contents (to various values of s') by the same solution, the value of C w is from definition 
independent of the water content (or e', equivalently). The errors of C w , represented by the dashed 
lines, depend on the standard error of determination of X s as a slope of the fitted straight line and on 
the standard error of determination of slope I of the X s vs. C s relation. 

The data points shown on the left panel of Figure 4 represent the values of C w calculated from 
Equation (11) with the use of the salinity index obtained from Equation (10), simulating field 
conditions (which may be named the "linear field model"). These values are obtained for each sample 
independently. The values of I were taken as in the previous linear model. The errors of C w , 
represented by the error bars, depend on the standard errors of all quantities used in Equation (10), as 
well as on the standard error of the parameter I. As can be seen, these errors are much greater than in 
the previous case, where X s from regression were used. Because of the form of Equation (10), as 
already mentioned, the errors may approach unreasonably high values for dry soils, which may be seen 
on the plots. Because of this, it was not possible to obtain reliable values of C w for the two samples 
with the least water content for each soil and moistening solutions; therefore, these values have not 
been presented in Figure 4. 

The values of C w obtained from the linear field model are roughly independent of the water content 
and agree with the values calculated from X s obtained from regression. A notable exception to this is 
soil no. 589 moistened with solution no. 4 (the most concentrated KC1 solution), as seen on the left 
panel of Figure 4. However, the values of C w calculated from X s obtained from regression are still 
encompassed within the error bars of the field model values. 

3.6. Quadratic Model of Salinity Index 

On the left panel of Figure 3, straight lines were fitted to approximate C b vs. e' relation for all 
measurement series. However, some possible quadratic effects may be lost. In order to assess the 
possible impact of quadratic terms on the C w values, a quadratic model of the salinity index has 
been examined. 

Even when the dependence C b vs. s' is not linear, it is still possible to use the salinity index defined 
as in Equation (6). With the quadratic relation: 

C b = a{s'f + be' + c (12) 
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which was fitted into the (e' , C b ) data points for all measurement series (not presented on figures), the 
salinity index calculated as a derivative takes the form: 

X s = 2ae'+b (13) 

and is obviously dependent on the real part of dielectric permittivity, and thus on the water content. To 
calculate electrical conductivities of soil water with the use of Equation (7), it is necessary to obtain 
slopes I of the X s vs. C s relations, for each value of e' . Fortunately, it transpired that the I vs. e' is 
linear (R 2 = 0.9988 for the worst fit for soil no. 568), which enables calculating C w for arbitrary 
values of e' . The values of electrical conductivity of soil water obtained with the use of the quadratic 
model of the salinity index are presented on the right panel of Figure 4. For comparison purposes, 
straight lines representing C w values calculated from X s obtained from regression from the linear 
model are also presented on the plots. One may notice that the values of C w either remain constant 
with s' (soil no. 605, 589, two series of soil no. 568) or decrease with the increase in s' . Also, for most 
measurement series, the values of C w obtained by all three approaches are comparable. 



Table 3. Conductivities of water content C w calculated from linear, linear field and 
quadratic salinity index models, with respective errors, for four moistening solutions, for 
samples of loess with similar values of the real part of the dielectric permittivity. 



No. of 
solution 


Re(e*) 


C w (linear model) 
(mS m" 1 ) 


Relative 
error (%) 


C w (linear field 
model) (mS m _1 ) 


Relative 
error (%) 


C w (quadratic 
model) (mS m _1 ) 


Relative 
error (%) 


1 


13.6 


263 


6.9 


278 


13.5 


255 


23.5 


2 


13.6 


758 


5.8 


801 


14.7 


756 


12.8 


3 


13.4 


1,424 


6.5 


1,434 


12.4 


1,441 


19.7 


4 


13.1 


1,919 


6.7 


1,796 


15.1 


1,921 


24.7 



The errors of C w calculated from the quadratic model, which depend mostly on the errors of the 
salinity index calculated from Equation (13), are quite high — with few exceptions, comparable to or 
even higher than for the linear field model, as can be seen in Figure 4 and Table 3. Table 3 shows a 
comparison of relative errors of soil water conductivity C w , determined from the salinity index 
obtained from: (i) regression from the linear model, (ii) calculated from Equation (10) (linear field 
model) and (iii) calculated from the quadratic model (Equation (13)). The presented errors pertain to 
C w obtained for loess for all four applied moistening solutions, for samples with similar values of the 
real part of dielectric permittivity. For other tested soils, these errors behave similarly. This suggests 
that the inclusion of the quadratic term in the C b vs. e' relation does not improve the C w estimation. 

Generally, only the linear field model may exhibit greater errors for C w than the quadratic one; 
however, not for all cases, as can be seen in Figure 4. Furthermore, because the salinity index in the 
quadratic model depends on water content, it would be very difficult to adapt it for field use. Because 
of this, the quadratic model does not present any significant improvement over the linear models for 
the tested soils. 
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3. 7. Conductivity of the Moistening Solution 

To further test the linear salinity index model, introduced by Malicki and Walczak in [7], one may 
calculate the conductivity of the moistening solution C s from Equation (9), by subtracting the residual 
conductivity C r (taken as conductivity of soil water for samples wetted with distilled water) from the 
conductivity of soil water C w : 

Cs ~ C w — C r , C r — C w (distilled water) (14) 

The results can be then compared with C s assumed from the molarity of the solutions (Table 2). The 
values corresponding to distilled water are equal to zero from definition. The results, calculated for the 
salinity index obtained as a slope of C b vs. s' relation (linear model), are presented in Figure 5. The 
values of R 2 and the standard errors cr presented in the figure describe how well the assumed values of 
C s explain the calculated quantities. As can be seen from the values, the agreement is very good. The 
worst case is soil no. 568, which exhibits both high errors in the calculated C s , as shown by the error 
bars, and the worst fit for the assumed values. The best case is soil no. 605, which has 95% sand 
fraction. This occurrence, along with the errors of C w presented in Figure 4, shows that the salinity 
index method is generally the most reliable for soils with the highest sand content. Because the 
calculated electrical conductivities of the moistening solutions have been determined from linear 
model, the values corresponding to distilled water are equal to zero from definition 

Figure 5. Conductivities of the moistening solutions calculated from the linear salinity 
index model vs. assumed conductivities for all tested soils. The straight line has a slope 
equal to 1 and corresponds to the perfect agreement. 



2000 



1 1500 ■ 



X 1000 • 



500 ■ 




500 1000 1500 2000 

C s assumed /mS m 1 



2000 



«g 1500 



1000 



500 




500 1000 1500 2000 
C s assumed /mS m 1 



2000 



<» 1500 




500 1000 1500 2000 
C s assumed /mS m 1 





2000 -1 






E 




1/1 


1500 ■ 




M 

X 


1000 ■ 


UIO. 




u 


500 ■ 




0 500 1000 1500 2000 
C s assumed /mS m 1 





2000 -J 






E 




i/i 


1500 ■ 


_E 




VI 




X 


1000 ■ 


E 




o 






500 ■ 


u 






0 ■ 




0 500 1000 1500 2000 

C s assumed /mS m 1 



♦ 



568 

589 

601 

605 

Loess 

V = « 



0.876, a=253.4 
0.972, 0=133.9 
0.995, o=61.31 



R 2 = 0.996, 0=51.16 
R 2 = 0.989, 0=85.15 



All soils: ^ = 0.971,0=85.84 



Sensors 2012, 12 



10904 



4. Conclusions 

The present study deals with the measurement of soil salinity using the FDR technique. It enables 
directly obtaining real and imaginary parts of complex dielectric permittivity. The obtained 
Equation (11) enables assessment of the soil pore water salinity from FDR measurements of soil bulk 
electrical conductivity and soil dielectric permittivity. Some of the parameters of this equation depend 
on soil properties, mostly clay content and specific surface area. Because of the coarse texture of the 
soil samples, those parameters have been determined just for sandy soils, and much research is needed 
for loam to clayey soils. In order to remove the influence of other physical factors and processes on the 
measured quantities, determination of water content was performed from the real part of dielectric 
permittivity in the frequency range 380-440 MHz, while the soil salinity was determined from the 
imaginary part of the dielectric permittivity for lower frequencies. 

The salinity index method, introduced in [7] for soil salinity estimation with the use of the TDR 
measurements, may still be applicable when the measurements are performed using the FDR 
technique. However, adaptation of the salinity index method for field use would require a good deal of 
work, because of the necessity to test wider varieties of soils with the application of more moistening 
solutions, in order to increase the accuracy of determining the electrical conductivity of soil water. 

The linear salinity model by definition gives electrical conductivity of soil pore water independent 
of the real part of dielectric permittivity and thus of water content. This ignores the possible ion 
adsorption and desorption processes occurring on the surfaces of the soil solid particles. Even though 
the chosen soils were sandy, they differed in specific surface area values (as shown in Table 1). 
Therefore, the application of the quadratic salinity index model, in which the salinity index depends on 
water content, could possibly detect adsorption effects. This was not so for the tested soils. However, 
the methodology presented in this paper will be tested in the future with soil samples with larger 
specific surface area values (loam and clayey soils). 
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